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ABSTRACT 

The transfer of polarized radiation in magnetized and non-magnetized relativistic plasmas is an 
area of research with numerous flaws and gaps. The present paper is aimed at filling some gaps 
and eliminating the flaws. Starting from a Trubnikov's linear response tensor for a vacuum wave with 
k = cj/c in thermal plasma, the analytic expression for the dielectric tensor is found in the limit of high 
frequencies. The Faraday rotation and Faraday conversion measures are computed in their first orders 
in the ratio of the cyclotron frequency f2o to the observed frequency lo. The computed temperature 
dependencies of propagation effects bridge the known non-relativistic and ultra-relativistic limiting 
formulas. The fitting expressions are found for high temperatures, where the higher orders in Q^/uj 
cannot be neglected. The plasma eigenmodes are found to become linearly polarized at much larger 
temperatures than thought before. The results are applied to the diagnostics of the hot ISM, hot 
accretion flows, and jets. 

Subject headings: radiative transfer — polarization — magnetic fields 
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1. INTRODUCTION 

We learn much of our information about astrophys- 
ical objects by observing the light they emit. Obser- 
vations of the polarization properties of light can tell 
us the geometry of the emitter, strength of the mag- 
netic field, density of plasma, and temperature. The 
proper and correct theory of optical activity is essen- 
tial for making accurate predictions. While the low- 
temperature pro pagation characteristics of plasma are 
well-established ([Landau fc Lifshit"ill98Clf ). the theory of 
relativistic effects has not been fully studied. In this pa- 
per I discuss the propagation effects through a homoge- 
neous magnetized relativistic plasma. A non-magnetized 
case emerges as a limit of the magnetized case. The dis- 
cussion is divided into three separate topics. 

Two linear plasma propagati on effects a re Fara day ro - 
tation and Faraday conversion (jAzzam fc Bashara 1987f ) . 
Traditionally, these effects are considered in their lowest 
orders in the ratio (3 of the cyclotron frequency JIq to the 
circular frequency of light w, id est in a high-frequency 
approximation. The distribution of particles is taken to 
be thermal 



dN = 



n cxp(— 7/T) 
47rm2r2if2(T-i)' 



P 



(1) 



with the dimensionless temperature T in the units of par- 
ticle rest mass temperature mc^/kB- The Faraday rota- 
tion measure RM and conversion measure are known in 
a non-relativistic T <^ 1 and an ultra-relativistic T ^ 1 
limits (Melrose 1997c). I derive a surprisingly simple an- 
alytic expression for arbitrary temperature T. 

The smallness of P = i^o/uj, P <C 1 in the real systems 
led some authors ()Melroselll997a|) to conclude that the 
high-frequency approximation will always work. How- 
ever, there is a clear indication that it breaks down at 
high temperatures T ^ 1. It was claimed that the eigen- 
modes of plas ma are linearly polarized for high temper- 
atures T 3> 1 ()Melroselll997cD . because the second order 
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term ~ f3'^ becomes larger than the first order term (3 
due to the T dependence. The arbitrarily large T-factor 
may stand in front of higher order expansion terms in /3 of 
the relevant expressions. I find the generalized rotation 
measure as a function of [3 and T without expanding in /3 
and compare the results with the known high-frequency 
expressions. The high-T behavior of the plasma response 
is indeed significantly different. 

Plasma physics involves complicated calculatio ns. This 
led to a number of errors in the literature (jMelrosd 
fT997^ . some of which have still not been fixed. In the 
article I check all the limiting cases numerically and an- 
alytically and expound all the steps of derivations. Thus 
I correct the relevant errors and misinterpretations made 
by previous authors, hopefully not making new mistakes. 
The analytical and numerical results are obtained in 
Mathematica 6 sy stem. It has a n enormous potential 
in these problems (jMarichevll2008f l. 

The paper is organized as follows. The formalism of 
plasma response and calculations are described in fj2l 
Several applications to observations can be found in 
fJH I conclude in 21 with a short summary and future 
prospects. 

2. CALCULATIONS 
2.1. Geometry of the problem 
I assume the traditional geometry depicted on Figure[TJ 

• Euclidean basis (e^, e^, e'^), 

• magnetic field along the third axis B ~ (0, 0, B)^ , 

• a wave vector of the wave k = fc(sin0, O,cos0)'^ 
with an angle 6 between k and B. 

The basis is rotated from (e^, e^, e'^) to (e^, e^, e"^), so 
that the wave propagates along k = (0, 0, k)^ in the new 
basis. The transformation has the form 



= cos sin ( 



2 -2 
e = e , 



= sin6'+e^ cos( 



(2) 



which can be conveniently written as 

/ cos 9 sin 9 



1 
sin 6 cos 9 



Vectors and tensors then rotate according to 



(3) 



(4) 




Fig. 1. — Geometry of the problem. 



2.2. Linear plasma response 

The propagation of weak electromagnetic (EM) waves 
in a homogeneous magnetized plasma can be fully de- 
scribed by the response tensor a'"'. It expresses the lin- 
ear proportionality between the induced current density 
and the vector potential j^{uj) — a^'^A'^{uj). The spatial 
projection of such defined 4-D tensor a'^i, is equal to the 
3-D tensor aij defined by j = aijA. 

I consider Trubnikov's form of the response tensor 
(jTrubnikov " 1958-. 'Melrose"1997a') . I work in a low-density 
regime, where the plasma response is calculated for a vac- 
uum wave with |k| = iv /c. I take the tensor g^" from th e 
first-hand derivations ijTrubnikovl 119581 : lMelrosdll997af ). 
make the transformation ([4]), and take the 1-st and 2-nd 
components in both indices. Thus the projection onto 
the (e^,e^) plane in CGS units is 



cmK2{p) Jo 



^ ( cos^ 9 cos 51o^ + sin^ 9 ij cos 9 sin JIqC 
' —rj cos 9 sin floS^ cos57o^ 



(5) 
(6) 



Rf" = ^^^^^ (cos 6l(sin - f^oO, ^^(1 ^ cos 17o0) 

"0 



Ri' = '^^^^ (cos ^(sin r^oC - f^oO: '7(1 - cos ^loO) > 

"0 

and 



(7) 
(8) 



p - 2iuj^p 



sin^ I 
"o 



1/2 



(9) 



where t] is the sign of the charge, Knir) is the n-th Bessel 
function of the second kind^ . The quantity p is the di- 
mensionless inverse temperature. 



(10) 



where Tp the actual temperature of particles. The re- 
sponse of plasma is usually characterized by the dielectric 
tensor. Its projection onto the (e^,e^) plane is 



(11) 



The wave equation for transverse waves in terms of ^ 
is 

E2 



0, 



(12) 



where E\ and E2 are the componen ts of the electri c field 
along and and — fc^c^/w^ ()Swansonll2003( ). 

2.3. High frequency limit 

Let me first calculate the limiting expression for a^^, 
in the high-frequency limit Qq <ti w. 1 denote 

no 



/3 = 



(13) 



substitute the definitions (fT5|) into the expression ([5]), 
and expand the response tensor a'^i, in (3. I retain only up 
to the 2-nd order of the expansion, whi ch gives the con - 
ventional generalized Faraday rotation ()Melroselll997cf) . 
The first terms of the series of r, t^^, and R^Ry read 



(5r , Tg = — 2iap^ Sr" — 



4«2 



R'^R^ 



12 ' ' 
(14) 

^ _ f l~cos^9-a'^f]'^/2 al3ricos9 \ , . 

^.sin^.(gj). (16) 

iMelrosd (|1997cf ) used the approximation Tq — —2iap 
instead of the expansion (fT4|l and obtained the approxi- 
mate high-T expressions as his final answers. 

However, one can take the emergent integrals, if one 
considers the exact expansions (|14ll5ll6p . Three terms 
appear in the expanded expression for a'^i, : 



da 



K2{ro) 



da 



da 



K'i{ro)5r 



21 



R^'R^ 



Kairo) 



(17) 
(18) 
(19) 



The 2-nd term (|T8|) originates from the expansion of 
K2{r)/r^ in r^ 



K2{r) K2{ro) _ Sr^ K^{ro) 



(20) 



'0 ^ '0 

Note that the analogous expression in IMelrosd Ill997cl) has an 
extra factor Oq? in the component t^^ and the opposite sign of 
Rf^R" term by an error. The author has corrected his formulas in 
IMelrosd l l2008t) . 
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Integrals (|17I18I19P can be evaluated knowing that 
da 



„K^Wp^-2ipa) 



2ipa 



p- 



(21) 



da 



,^ K,i^p^~2ipa) 



^,^„+iX^^_ (22) 



2.4. Components in high-frequency limit 

I substitute the high-frequency expansions (|14I15I16P 
into the expression (fTTj) for the projection of the dielectric 
tensor e^i, with the projection of the response tensor a^i, 
(O and take the integrals (|17I18I19|) analytically. The 
components of the dielectric tensor pT|) in the lowest 
orders in ^Iq/uj are then 



a;2 \K2{p) V ^ 



LLp-p 



(23) 

4 f^l<£)C,^2|)+I2lil!f). ,24) 



if2(p) 

where the plasma frequency in CGS units is 
2 47rng2 



(25) 



(26) 



The results reproduce the non-relativistic limits for 
p — > +00 : 



£^-l-4fl + ^COs20|, 



2 . "^p^O 



e 2 = — £ 1 = — 



cosy, 



(27) 
(28) 
(29) 



where all Bessel functions of p approach unity 



. ap j: ., 

(iLandau fc Lifshitslll980l : lTrubniko\ll996l : ISwansonl2003l : 
lBellanl[2006l ). The corresponding relativistic limits p ^ 
of the same components are 

e\ = 1 - ^ ( — ( 1 + ^ cos^ Q 1 + 



cj2 V 2T V w 



cj2 V 2T V w 



t^,-? / 1 / fi^?\ „7f^^sin2' 



e 2 = -£"1 = -i?? 



i^l^o ln(r) 



cosy, 



(30) 
(31) 

(32) 



t£ 2T2 

cons istent with iMelrosd ()1997cD : iQuataert fc Gruzinovl 
(|2000f )^. The ultra-relativistic non-magnetized disper- 
sion relation then reads 

2 I 2,2 27rng2 

= — ^ + c fc = V c k 

2T mT 



(33) 



The non-diagonal term has a wrong sign in IMehosd a997d). 
^ Th e diagonal plasma response is 2 times larger in IMelrosd 
Ill997d ') by an error. 



according to t he relation (fT2l). The exp ression (|33p is 
consistent with ILandau fc LifshitsI ()198Clf ). chapter 32. 

The plasma propagation effects can usually be de- 
scribed in terms of only the difference of the diagonal 
components and the non-diagonal component of e^^. I 
define X to be a vector of T, 6, flo/oj. I introduce the 
multipliers /(X) and (7(X) to correct the expressions, 
when the high-frequency limit breaks. I write the differ- 
ence between the diagonal components with a multiplier 
/(X) as 

e\^e^ = fiX)^(^l^+6T)sin^e (34) 



V^2(r-i) 

and the non-diagonal component with a multiplier (?(X) 



as 



^^»oXo(T-i) 
uj3 K2{T-^) 



COS! 



(35) 



Both multipliers equal unity in the high-frequency limit 
/(X) — 5(X) = 1. Now we can turn to a more general 



2.5. Fitting formulas for higher temperatures 

The ultra-relativistic expressions pOI31l32p allow me 
to trace the T-factors in front of the first 3 expansion 
coefficients of the dielectric tensor in (3. The coefficient 
at 0^ is ~ T3/ln(T) times larger than at p. Thus at 
temperature T > 10 the 2-nd order becomes larger than 
the 1-st order for the ratio VIq/lo ^ 10"'^. This indicates 
that the expansion in (3 may become invalid at these 
plasma parameters'*. The multipliers /(X) and (?(X) are 
likely to be far from 1. I consider only the real parts of 
these multipliers, since the imaginary parts correspond 
to absorption. The contour plots of the numerically cal- 
culated /(X) and (7(X) for somewhat arbitrary 6 = tt/4 
are shown on Figure [5] and Figure [H respectively. 




-3.0 -2.5 

Log,„[!i„/a 



Fig. 2. — Multiplier /(X) for the difference of the diagonal 
components £-"^1 — e^2 for S = 7r/4. 



One cannot claim that th e diagonal mag netized terms become 
larger then the non-diagonal l lMelrosdll997dV 



Fig. 4. — Multiplier f{X) for the difference of the diagonal 
components e^i — e^2- Dashed line — fitting formula I I37I I. 
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Fig. 3. — Multiplier (/(X) for the non-diagonal component e^2 
for e = 7r/4. 

Let me define X to be the following combination of the 
parameters 



X ^TjV2sm9 i 103 — 

UJ 



(36) 



For the fiducial flo/uj = 10^'^, 6 = -k/A the parameter X 
is just temperature X = T. 

I first identify the boundaries, where the high- 
frequency limit is valid. Then I find a fit for the multi- 
pliers at higher X. The expression p4p for the difference 
e^i — is accurate within 10% for X < 0.1 if we set 
f{X) = 1. The expression (|55|) for £-^2 is accurate within 
10% for X < 30 if we set g{X) = 1. The accuracy de- 
pends on the parameter X rather than on the individual 
parameters T, D,q/lu, 9. The expression 



f{X) = 2.011 exp - 



X 



1.035 



4.7 



/ X 

0.011 cxp 

2.73 / ^ V 47.2 



(37) 



' cos [ y ) exp 



extends the applicability domain of the formula (p4|) up 
to X 200. Figure m shows the fit for f{X) in compari- 
son with the numerical results. The expression 

g(X) = l-0.111n(l + 0.035X) (38) 

extends up to X ^ 200 the domain of the formula ([35]) . 
Figure shows the fit for g{X) in comparison with the 
numerical results. 



2.6. Exact plasma response 

The expression for the response tensor ([5]) is written 
for a vacuum wave with |k|c = uj. In the real plasma, 
the wave is modified by the plasma response. A more 
;eneral self -consis tent resp onse tensor should be used 
Trubnikovl [l958l : iMelrosd Il9973 l. One needs to solve 
a dispersion relation similar to the relation (fT2|) to ob- 
tain the eigenmodes. Thus the eigenmodes and the re- 
sponse tensor should be computed self-consistently. One 
should not forget about the antihermitian and longitudi- 
nal components of the dielectric tensor e^i, that modify 
the dispersion relation. 




Fig. 5. — Multiplier g{X) fo r th e non-diagonal component £^2- 
Dashed line — fitting formula l|38| l. 

2.7. Eigenmodes 

The above calculation is applicable also to a non- 
magnetized plasma. Dispersion relation of EM waves in 
a non-magnetized plasma reads 



7 2 2 

k c 



•K2{T- 



(39) 



in a high-frequency approximation uj ^ Th e oppo- 
site limit oi kc <ti to was considered by 'Bergman ( 2001') . 

Now we turn to the magnetized case. llVIelrosc ( 1997d ) 
only considered the first terms of in the expansion of 
a^jy in (3 to get the eigenmodes. I do the next step: con- 
sider the full expression in /? in the low-density regime 
kc — LU, but consider only the hermitian part of a'^^ in 
computations. The ellipticity T = (e^i — £^2) : \£^2\ de- 
termines the type of eigenmodes. If |T| ^ 1, then the 
eigenmodes are linearly polarized unless 9 is close to 0. If 
|T| <C 1, then the eigenmodes are circularly polarized for 
9 far from 7r/2. Let mc consider the fiducial model with 
flo/u! = 10~3 and 9 = 7r/4. Figure [6] shows the ratio T 
calculated in a high-frequency approximation (see § 12. 3|) 
(dashed line) and in a general low-density approximation 
(see ? 12. 5p (solid line). The high-frequency approxima- 
tion produces the linear eigenmodes already at T > 10 
consistently with , Melrose (i1997g ). However, the general 
low-density limit produces the eigenmodes with T ~ 1 
up to very high temperatures T ~ 50. Unexpectedly, 
the sign of the diagonal difference (e^i 
about T « 25. 

3. APPLICATIONS 



£ 2) changes at 




in hot magnetized plasma 



1.0 1.6 2.7 4.4 7.3 12. 20. 33. 54. 90. 148 244 403 
Temperature T 

Fig. 6.— EUipticity T = (e^i - e^j) : I^IjI of eigenmodes. The 
absolute value of the ratio T much above unity — linear eigen- 
modes, much below unity — circular eigenmodes. Solid line — 
this paper, dashed line — previous calculations. 

The calculated transrelativistic propagation effects 
have far-reaching consequences in many topics of astron- 
omy. Let me concentrate on four applications: propa- 
gation delay, Faraday rotation measure of light from the 
Galactic Center (GC), circularly polarized light from the 
GC, diagnostics of jets. 

3.1. Dispersion measure 

Propa gation delay is an important effect in pulsar dis- 
persion ( Phillips fc WolszczanI |1992[ ). The relativistic 
part of this delay can be obtained from the dispersion 
relation (|39p . I retain only the first-orde r correction in T, 
since T ^ 1 in the interstellar medium (|Cox k, Reynolda 
[1981 . Since Ki{T-^) / K2{T-^) « 1 - 2,T/2 at low T, 
the non-relativistic Dispersion Measure (DM) should be 
modified as 



DAL, 



-T 



(40) 



This shows that the gas density is slightly underesti- 
mated, if the non-relativistic formulas are used^. How- 
ever, the relativistic correction to the DM is small and 
can be neglected in most practical cases when T ^ 1. 
The effects in magnetized plasma are also relevant for 
pulsars. 

3.2. Magnetized radiative transfer 
3.2.1. General formulae 

Relativistic plasmas exhibit a generalized Faraday ro- 
tation for a general orie ntation of the magnetic field 
(jAzzam fc Basharalll987() . One can decompose it into 
two effects: Faraday rotation and Faraday conver- 
sion. The former operates alone at = 0, tt, the lat- 
ter operates alone at 9 = 7r/2, and both should be 
considered together for the intermediate angles. The 
transfer equations (Mueller calculus) for the Stokes 
parameters I, Q, U, V were devised to treat to- 
gether the propagat ion ef fects, emission, and absorptio n 
(jAzzam fc basharal [l987l : iMehose fc McPhcd ranI Il991f ). 
Good approximat ions for emission and absorption have 
been long known (iTrubn ikov 1958; Rybicki fc Lightmai 
[19671 : iMelrose fc McPhedranlil991t , Wolfe fc Meha.200 ^ 
Now one can combine them with the proper approxima- 
tions of the propagation effects given by 






-pv pu 

Pv -pQ 

-pu PQ 




(41) 



Pv 



-te^2, PQ = -^{e\-e\), pu ^ 0, (42) 



and do the radiative transfer calculations. Here e'^^ 
stands for the Hermitean part given by the relations 
(|37I38P with the real multiphers f{X) and g{X). One 
of the most interesting objects for such calculations is 
our Galactic Center Sgr A*. 

The transfer equations were recently solved for a sim- 
ple tim e- independent dyna mical model of the GC ac- 
cretion (jHuang et al.ll2008l ). The authors treat the or- 
dinary and extraordinary modes as linearly polarized. 
They assume these eigenmodes constitute a basis, where 
cither U or Q components of emissivity and propaga- 
tion coefficients vanish. Actually, U components vanish 
{pu = 0) already in the basis (e^,e^), since the pro- 
jection of the magnetic field onto (e ^,e^) is parallel to 
fsee IMelrose fc McPhedranI ^91) p. 184). As I have 
shown in the § 12. 7( plasma modes are far from being 
linearly pola rized at temperatur es T < 10 estimated 
for the GC (jSharma et al.ll2007f) . Thus, the propaga- 
tion coefficients should be taken from equations (|34|) and 
(|35|) . The Faraday conversion coefficient pq cannot be 
defined via emissivities and Faraday rotation coefficient 
Pv as in lHuang et al.l ()2008[ ). The Faraday rotation mea- 
sure was calculated from a s imula ted accretion profile in 
ISharma. Ouataert fc Stoni ()2007D . However, the paper 
considered only the Faraday rotation and did not carry 
out the self-consistent treatment of propagation. It is 
impossible to disentangle the effects of Faraday rotation 
and Faraday conversion in a relativistic plasma. 



3.2.2. Faraday rotation 

The crucial part of any radiative transfer is the 
proper transfer coefficients. It allows one to esti- 
mate the electron density near the accreting object 
([Ouataert fc Gruzino v 2000; .Shcherba kov 200g)- Sev- 
eral formulas were suggested for the temperature de- 
pendence of the component s^2 responsible for Fara- 
day rotation. These formulas were yet given for the 
high-frequency approximation (see § 12. Sp . Let me com- 
pare them with the exact temperature dependence (I25p 
J = Kq{T-^)/K2{T-'^) and its limits. The limits are 
J ^ 1 as T -> and J ^ ln(T)/(2r2) as T ^ +oo. The 
results of this comparison arc shown on Figure [7| 



0.16 



■D 0.016 



0.005 




0.16 0.4 1. 2. 5. 
Temperature T in (m c^/kg) 



5 The formula in IPhillips fc WolszczanI llT99l) has no refer- 
ences/checks and is not correct. 



Fig. 7. — Temperature dependence of the Faraday rotation mea- 
sure. 



6 



iBallantvne. Ozel fc Psaltii (|2007[ )'^ divided the 
thermal distribution into ultra-relativistic and non- 
relativistic parts as marked by the electron energy 
7crit = 10. They sum the contributions of both species 
with calculated densities. To make a plot, I take their 
effective temperature Q of plasma above 7crit to be 
just tempe rature 8 = T and not th e average kinetic 
energy as IBallantvne. Ozel fc Psaltii (|2007D suggest. 
This brings O to lower values and decreases the rotation 
measure. Even with this decrease the rotation measure 
is severely overestimated at T ~ 1. The convergence to 
the relativ istic limit is not a chieved even at T '--^ 30. 
The paper [Huang et aO (|2008l ) found the simpler fitting 
formula that reproduces the limits. Their expression is 
quite accurate.^ 

3.2.3. Faraday conversion 

The increase in the circular polarization of Sgr A* at 
frequency ITHz is predicted by[Huang ct al. (2008). The 
phase of Faraday conversion approaches unity and the 
destructive interference does not occur at this frequency. 
The result seems to be qualitatively correct regardless of 
the expression for the conversion measure, but the proper 
expressions and ([55|l should be used for quantitative 
predictions. 

3.2.4. Jets 

The better treatment of propagation effects may also 
play a role in observations of jets. As we saw in § 12. 5( 
the propagation effects in thermal plasma cannot be de- 
scribed in the lowest orders in ^o/^j if the temperature 
T is sufficiently high. Power-law distribution of elec- 
trons can have a quite high effective tempera ture. Thus 
the high-frequency li mit (Sazonoy 1969; Jones fc O'Delll 
Il977t lMelrosel[l997b[ ) may not approximate well the her- 
mitian part o f the response tensor. Careful analysis of jet 
observations (jBeckert fc Falckdl2002l: iWardle et aO 19981) 
may be needed. It should be based at least on the ex- 
pressions for e^u in a general low-density regime. 

^ The paper IBallantyne. Oz el &: PsaltisI l|2007l ') has likely con- 
fused the 3-D projection of the 4-D response tensor in = o^^A^ 
Ipvlehose 1997c ) with the 3-D response tensor j = OijA that has 
the opposite sign. 

" Temperature" 7c in IHuang et al.l l|2008 ) should be redefined 
as 7c = 1 -f T, otherwise the lower limit is not reproduced. 



4. DISCUSSION & CONCLUSION 

This paper presents several new calculations and 
amends the previous calculations of propagation effects 
in uniform magnetized plasma with thermal particle dis- 
tribution equation The expression ([5]) for the cor- 
rect response tensor is given in a high-frequency ap- 
proximation. The exact temperature dependence (|2.4p 
and (|25p is found in first orders in flo/cu in addition to 
the known highly-relativistic and non-relativistic results. 
The higher order terms may be important for relativis- 
tic plasmas in jets and hot accretion flows. The fltting 
expressions (1571) and ([55)1 are found for the dielectric ten- 
sor components ([34]) and ([35]) at relatively high temper- 
atures. 

The results of numerical computations are given only 
when the corresponding analytical formulas are found. 
One can always compute the needed coefficients numeri- 
cally for every particular frequency lu, plasma frequency 
Wp, cyclotron frequency fio, and distribution of electrons. 
However, the analytic formulas offer a simpler and faster 
way of dealing with the radiative transfer for a non- 
specialist. The eigenmodes were not considered in much 
detail, since radiative transfer problems do not require 
a knowledge of eigenmodes. However the knowledge of 
eigenmodes is needed to compute the self-consistent re- 
sponse tensor (see § 12. 6p . 

The response tensor in the form (O can be expanded 
in Qq/lu and ujp/uj. This expansion is of mathematical 
interest and wiU be presented in a subsequent paper as 
well as the expressions for a power-law electron distribu- 
tion. Propagation through non-magnetized plasmas will 
also be considered separately. 



The author is grateful to Ramesh Narayan for fruitful 
discussions and Diego Munoz for pointing out relevant 
references. I thank the anonymous referee for helpful 
suggestions that improved the paper. 
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